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ABSTRACT 

Many extrasolar planets orbit closely to their parent star. Their existence raises the fundamental 
problem of loss and gain in their mass. For exoplanet HD209458b, reports on an unusually extended 
hydrogen corona and a hot layer in the lower atmosphere seem to support the scenario of atmospheric 
inflation by the strong stellar irradiation. However, difflculties in reconciling evaporation models 
with observations call for a reassessment of the problem. Here, we use HST archive data to report 
a new absorption rate of r^ 8.9% ± 2.1% by atomic hydrogen during the HD209458b transit, and 
show that no sign of evaporation could be detected for the exoplanet. We also report evidence of 
time variability in the HD209458 Lyman-a flux, a variability that was not accounted for in previous 
studies, which corrupted their diagnostics. Mass loss rates thus far proposed in the literature in the 
range 5 x 10^° — 10^^ gs"'^ must induce a spectral signature in the Lyman-a line profile of HD209458 
that cannot be found in the present analysis. Either an unknown compensation effect is hiding the 
expected spectral feature or else the mass loss rate of neutrals from HD209458 is modest. 
Subject headings: Stars: individual (HD209458) — Stars: planetary systems — Ultraviolet: stars — 
Line: profiles — Techniques: spectroscopic 



1. INTRODUCTION 

Of all the planets discovered outside the realm of our 
solar system, some of the most dramatic new classes of 
objects are those in which the planet is a gas giant or- 
biting at merely a few stellar radii (^ 0.02AU) from its 
parent star. These close-in extrasolar planets are Jupiter- 
like giants that are exposed to strong fluxes, magnetic 
fields, and plasma winds — a very harsh and active stellar 
environment. Because of their stars' proximity, gravity, 
through tidal effects, distorts the shape of their atmo- 
sphere while the continuous extreme ultraviolet (UV) 
energy deposition inflates it (Lammer et al., (2003); 
Baraffe et al., (2004); Lecavelier et al., (2004); Yelle 
(2004); Jaritz et al., (2005); Tian et al., (2005); Munoz 
(2007); Villaver and Livio (2007)). Unfortunately, little 
is known about those regions that separate extrasolar gi- 
ant planets from their stars, particularly the immediate 
environment of the planet. 

One of the most extensively studied extrasolar systems 
is HD209458. For reference, HD209458b was flrst discov- 
ered transiting its parent star and covering 1.5% of its 
disk (Charbonneau et al., (2000); Henry et al., (2000)). 
Some of the flrst attempts to learn about the immediate 
environment of the planet were Lyman-a (121.6 nm) ob- 
servations of the system using the Space Telescope Imag- 
ing Spectrometer (STIS) onboard the Hubble Space Tele- 
scope (HST). A flrst program of observation, obtained 
with the STIS/G140M grating and the 52"x0.1" slit, was 
implemented in 2001 during HD209458 planetary transit 
but no conclusions were reported (see Table 1). Soon 
after, a second program visited the target during three 
transits (Vidal-Madjar et al., (2003)). An initial analysis 
of this data set concluded that a huge cloud of hydrogen 
is covering 15% ±4% of the stellar disk (Vidal-Madjar et 
al., (2003)); it also claimed that spectral absorption dur- 
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ing transit is deeper on the blue side of the stellar line. 
Accordingly, the hydrogen cloud was required to extend 
beyond the planetary Roche limit where an intense es- 
cape of ~ 10^*^ gs^^ of hydrogen is a priori operating. 
These results, along with other far UV low-resolution 
observations of heavy constituents, led to the conclu- 
sion that the upper atmosphere of HD209458b should 
be in a hydrodynamic blow-off state (Vidal-Madjar et 
al., (2003); Vidal-Madjar et al. (2004)). 

Numerous studies then followed on different mecha- 
nisms for hydrogen loss from hot exoplanets closely or- 
biting their stars (Lammer et al., (2003); Baraffe et 
al., (2004); Jaritz et al., (2005); Lecavelier et al., 
(2004); Yelle (2004, 2006); Tian et al., (2005); Munoz 
(2007)). As noted by Munoz (2007), all loss rates 
thus far proposed by theoretical models in the range 
5 X 10^° — 10^^ gs~^ exceed the lower limit provided by 
Vidal-Madjar et al., (2003). Unfortunately, most studies 
neglected to quantitatively translate their loss rate to a 
spectral absorption in the Lyman-a line proflle that could 
be tested with the HST observation. Independently, 
pointing out that the observed mass function distribution 
of extrasolar giant planets (EGPs) follows a trend M~^ 
for mass range ^ 0.2 — 5 Mj, where Mj is the Jovian mass, 
Hubbard et al., (2007) derived the same mass function 
distribution for highly irradiated EGPs orbiting at dis- 
tances smaller than ~ 0.07 AU. Accordingly, Hubbard et 
al., (2007) rejected substantial mass loss during EGPs' 
migration to smaller distances from their star, unless the 
loss mechanism is compensated by an unknown process. 
When combined with the unusual scales derived for the 
hydrogen extent and escape, all these studies then call for 
a careful reassessment of the HST Lyman-a observations 
thus far obtained on HD209458, at least in order to pro- 
vide validated constraints on theoretical models. 

2. OBSERVATIONS AND DATA ANALYSIS 
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In the following, we report a new analysis of archive 
data obtained during the two HST/STIS programs de- 
scribed above. In total, we have four visits of the tar- 
get corresponding to three exposures of roughly ^ 2000 s 
duration each, resulting in 12 exposures of the systems 
around the transit period (Table 1). All observations 
were obtained in the time-tag mode, a technique that 
keeps track every 125 x 10~^ s of photon events during 
each exposure. The question, then, is: why is this mode 
of observation important in the present case? First, we 
stress that the transit effect is a weak variation of the 
stellar signal. As such, its trend is best represented by 
a dense time series. Second, chromospheric and coro- 
nal variabilities of the star are unknown in the Lyman- 
a spectral window considered here, and this may seri- 
ously corrupt any diagnostic. 

To properly handle time tagged data, we partition each 
time-tag exposure into a set of shorter sub-exposures, 
taking into account the heliocentric and barycentric time 
correction procedure in which ephemerides are retrieved 
from the archives before the IRAF/STSDAS "odelay- 
time" procedure is apphed (Brown et al., (2002)). After 
several trials, we found that 300 s sampling of the time- 
tag data is a good compromise between acceptable signal- 
to-noise ratio (S/N) and time coverage. Next, each sub- 
exposure is calibrated through the STIS pipeline. The 
emitting source (the star plus sky background) is pre- 
sumed extended, an option that allows efficient control 
of the subtraction of the sky background contamination. 
Time is then converted to fixed orbital phases measured 
from the transit central time (TCT), itself carefully taken 
from the most recent and accurate determination of the 
HD209458 system parameters (Ballester et al., (2007); 
Knutson et al., (2007)). Sub-spectra of identical phase 
positions are accumulated from the initial twelve expo- 
sures, resulting in a unique 53 bins time series of the 
system versus the orbital phase angle (Fig ??). Unex- 
pectedly, three gaps, lasting respectively 312 s, 404 s, and 
406 s, appear in the time series, for which no observa- 
tion is available. Because the three gaps are narrow and 
well-separated from each other, we determined that fill- 
ing them has a negligible effect on our final conclusions 
(Schneider (2001)). Errors due to photon counting have 
been propagated, taking into account the correlation be- 
tween the different phase positions relative to the initial 
sampling of sub-exposures time over the full observing 
time period. 

We next define the wavelength domain of contamina- 
tion by the sky background, including both the Earth's 
geocorona and the interplanetary medium emissions. 
The difficulty comes from the uncertainty about the geo- 
corona's strength when estimated from a detector sec- 
tor, along the STIS slit, different from the one where 
the stellar signal was recorded. First, we subtracted the 
dark noise of the detector following (Vidal-Madjar et al., 
(2003) and Ballester et al., (2007)) and then compared 
the sky background signal from different sectors along 
the slit and for different conditions of observation. Our 
conclusion is that the STIS MAMMA detector has an in- 
herent non-uniformity corresponding to an incompress- 
ible uncertainty of 5% on extended sources. Coinciden- 
tally, this uncertainty is comparable to photon statistical 
errors. To ensure that such error will not corrupt the 
stellar signal per wavelength pixel at the 1% level, we 
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Fig. 1. — Contours of equal flux for the HD209458 time series 
vs wavelength (nm) and time from transit (or planetary orbital 
phase angle). Top axis shows velocities in the stellar rest frame 
(~ lOkms"^) from heliocentric system). Horizontal dark bands 
are time gaps of ~ 300-400 s for which no observation is available. 
These bands are linearly interpolated using nearby spectra for light 
curves analysis and in/out transit spectra comparison. Following 
the time evolution of the signal (from bottom to top) , we observe 
the transit absorption as a slight dimming of the stellar flux start- 
ing ~ 5000 s before TCT. The wide, vertical, light purple band 
corresponds to the spectral window of extinction of the stellar sig- 
nal by the interstellar gas along the line of sight (Wood et al., 
(2005) 

deduce that a void window [121.541, 121.584] nm should 
be disregarded in any spectral analysis that requires high 
accuracy, such as for a transit event or short-term stellar 
variability. 

3. LIGHT CURVE TREND AND TIME VARIABILITY 

Next, to obtain the trend of the planetary transit 
in Lyman-a , we integrated the time series spectra in 
the range [121.483,121.536] nm on the blue wing and 
[121.589,121.643] nm on the red wing. These ranges 
were selected so that the stellar signal per wavelength 
pixel remains above the statistical noise (> Ict). The 
resulting light curve is noisy, but a trend is apparent and 
can be efficiently extracted (Fig. 2a). To dampen the 
signal noise while keeping a clear trend, the best compro- 
mise is to gather data by eight phase bins for a new bin 
of 2400 s. Accumulating the signal from three new bins 
(3 X 2400 s) inside transit, we derive ~ 8.9 ± 2.1% drop- 
off of the stellar Lyman-a intensity during the planetary 
transit (Fig. 2a). Our absorption rate of ~ 8.9% ± 2.1% 
is much lower and accurate than reported in a previ- 
ous study (Vidal-Madjar et al., (2003)), yet a marginal 
agreement could be found between our maximum rate 
(11%) and their bottom value. If this obscuration is con- 
verted directly to a planetary occulting disk, then one 
would obtain a hydrogen cloud of ~ 2.47 ± 0.30 Rp ra- 
dius, much smaller than the Roche lobe limit of ~ 4.08 
Rp (Gu et al., (2003)), where Rp = 1.32 Rj is the most 
recent estimate of the radius of IID209458b (Ballester et 
al., (2007); Knutson et al., (2007)), and Rj is the Jo- 
vian one. Now, to capture the trend of the transit curve, 
we used a sophisticated 2D model of planetary transit 
at Lyman-a that accurately accounts for the atmospheric 
radial structure of the planet (Yelle (2004); Ben-Jaffel et 
al., (2007)) and properly estimates the atmospheric ob- 
scuration versus wavelength, including extinction by the 
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TABLE 1 

HST/STIS Data Set on HD209458 Used in This Study. All observations were obtained with the G140M grating and the 

52"x0.1" long slit. Transit central time (TCT) is defined by 2,452,826.628521 HJD (Ballester et al., (2007); Knutson et 

AL., (2007)).' 



Dataset name 


Program ID 


Start time - TCT (s' 


Duration (s 


End time - TC 


O4ZEA4010 


7508 


-7980.28 


2600. 


-5380.12 


O4ZEA4020 


7508 


-2202.29 


2600. 


397.87 


O6E201010 


9064 


-8075.11 


1780. 


-6295.02 


O6E201020 


9064 


-2246.12 


2100. 


-145.96 


O6E201030 


9064 


3531.89 


2100. 


5631.98 


O6E202010 


9064 


-10167.67 


1780. 


-8387.50 


O6E202020 


9064 


-4750.68 


2100. 


-2650.52 


O6E202030 


9064 


1025.31 


2100. 


3125.46 


O6E203010 


9064 


-11258.45 


1780. 


-9478.27 


O6E203020 


9064 


-5876.47 


2100. 


-3776.29 


O6E203030 


9064 


-102.45 


2100. 


1997.73 
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Fig. 2. — (a) Light curve (LC) obtained from our time scries (as 
shown in Fig. 1) by integration of the signal in the spectral windows 
[121.483, 121.536] nm and [121.589, 121.643] nm. Histogram and 
related errors (plotted every two bins for clarity) show the LC 
with a time bin of 300 s, while filled circles with attached error bars 
represent the LC rebinned to a larger timescale of 8 X 300s = 2400s. 
The solid curve is our best least square fit to the rebinned LC. A 
total obscuration of ~ 8.9% it 2.1% is derived during the planetary 
transit, (b) Ratio of LC (histogram) to best model fit (solid line) 
is shown with related statistical error bars. This ratio cancels the 
transit trend shown in (a). The resulting signal is a good indicator 
of the variability of the HD209458 Lyman-a intensity vs time. 

interstellar gas intervening along the line of sight (Wood 
et al., (2005)). Our best least square fit is shown in Fig- 
ure 2a. For our purpose of time analysis, we remark that 
a functional fit could also be a good model to obtain the 
light curve's trend. 

We can now determine the stellar signal time evolu- 
tion after we cancel the transit trend using our best fit 
to the observed light curve (Fig. 2a). The resulting ra- 
tio shows a variable behavior with an average amplitude 
^ 8.6 ±5.6% of the stellar integrated intensity (Fig. 2b). 
Using the Durbin- Watson statistical test (Durbin and 
Watson (1951)), we found no apparent serial correlation 
at the 1% confidence level in the corrected signal — a sig- 
nal that also shows no evident periodicity. HD209458 
was previously suspected to have a relatively moderate 



chromospheric activity from Call H and K lines that were 
recorded over full orbits of the system (Shkolnik et al., 
(2005)). Our finding of a time variation of 8.6% on aver- 
age in the stellar Lyman-a signal, with peaks that may 
reach ^ 20% (> 3cr ), seems to support a relatively active 
corona of the star, presumably up to the planet's orbit. 
Such activity could be of common origin (flaring, non- 
uniformity of the stellar disc during transit, etc.) and/or 
related to an enhancement of magnetic activity on the 
star-planet line (Zarka (2007)). Also, one can speculate 
about the hydrogen cloud topology around the planet 
and its evolution with time. To that end, comparative 
studies with interacting binary stars may be useful in 
clarifying the different regimes of interaction between an 
exoplanet and its host star (Shore et al., (1994)). Un- 
fortunately, the FUV observations thus far obtained do 
not cover a full orbit of the planet, thereby making it 
difficult to predict the exact configuration of the star- 
planet system. In any case, we believe that the unusual 
15% obscuration previously reported (Vidal-Madjar et 
al., (2003)) was corrupted by this unaccounted-for vari- 
able component in the star-planet system signal. Here we 
can extract it because we are able to sample the transit 
period by a dense time series using the information gath- 
ered from the time tag mode of HST/STIS and ~ 25% 
more observation time from the archives. 

4. PLANETARY MASS LOSS OR FLUX VARIABILITY? 

In the following, we compare the in/out of transit stel- 
lar line profiles. The impetus of this study is the need 
to determine the relevance of a blucshiftcd absorption in 
the stellar line profile that may occur during transit, as 
claimed in earlier studies (Vidal-lVIadjar et al., (2003)). 
On the one hand, we derive an average unperturbed pro- 
file of the HD209458 Lyman-a emission line by merging 
all sub-spectra of the time series that we correct for the 
transit trend with the best fit shown in Figure 2a. The 
resulting profile is a good reference that best represents 
the out-transit stellar line and for which time variability 
has been reduced to the 1% signal level (Fig. 3a). On 
the other hand, the in-transit line profile, when corrected 
for the ~ 8.9 % drop-off during transit, properly recov- 
ers the unperturbed line profile (Fig. 3a), leaving no real 
possibility of extra absorption as claimed in prior studies 
(Vidal-JVIadjar et al., (2003)). 

To further investigate how time variability of the 
HD209458 Lyman-a emission line corrupts the diagnos- 
tic as it pertains to extra absorption or emission features 
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Fig. 3. — Comparison between Lyman-aline profiles in and out 
of transit period. The sky background spectral window is indicated 
by two dashed vertical lines, (a) The in-transit line profile (solid 
thin line) is accumulated for the time period starting ^ 3900 s 
before TCT and ending ~ 3900 s after it. To correct for the ~ 8.9 
% obscuration derived in this study, the corresponding intensity is 
scaled by 1.098. The resulting line profile (dotted curve) properly 
recovers the unperturbed line profile (histogram), (b) The first 
in-transit line profile, Bl (thin solid line), was accumulated over 
the time period starting ~ 4000 s before TCT and ending ~ 600 s 
after it. The second in-transit line profile, B2 (dotted line), was 
accumulated over the time period starting ~ 1800s before TCT 
and ending --^ 3900s after it. 

that may appear in the stellar line during transit, we se- 
lected two phase windows inside the transit period for 
which we compared line profiles to the unperturbed stel- 
lar line. As shown in Figure 3b, a direct comparison 
would indicate that line peaks are equally absorbed for 
line profile Bl, while for line profile B2, the red peak is 
the most absorbed. On the basis of line profile B2, the 
diagnostic would be just the opposite of that of (Vidal- 
Madjar et al., (2003)), leading to escaping hydrogen 
toward the star, while for line profile Bl, the diagnostic 
would be no H escape. The problem is that these inter- 
pretations of preferred blueshifted or redshifted absorp- 
tion do not account for the relatively strong modulation 
of the stellar signal evidenced in this study. Therefore, 
any claim of a preferred absorption during transit, either 
blue or redshifted, is not realistic, particularly at this rel- 
atively modest level of the signal to noise. It follows that 
the blueshifted absorption, advanced in previous studies 
(Vidal-Madjar et al., (2003); Lecavelier et al., (2004)) 
as a signature of atmospheric evaporation in a cometary- 



like tail of HD209458b, has, unfortunately, no foundation 
in the HST/STIS data set as it was only the effect of the 
stellar signal variability with time that corrupted the di- 
agnostic. 

5. CONCLUSION 

We use HST archive observations of the Lyman- 
a emission of HD209458 to report an absorption rate of 
^ 8.9% ± 2.1% by atomic hydrogen during the transit 
of the planetary companion. If the planet is sketched 
as a compact blocking body, our analysis requires an H 
cloud effective extent that does not exceed ^2.5 Rp — 
a size that falls short of the Roche limit r^ 4.08 Rp of 
IID209458b. In addition, time variability of the stel- 
lar flux is evidenced, but no sign of extra or Doppler- 
shifted absorption could be detected during transit. This 
absence of extra absorption during transit and the rel- 
atively small size of the effective area of the hydrogen 
cloud around the exoplanet make it difficult to conceive 
of significant atmospheric evaporation from the planet. 
Of course, we cannot rule out that a complex atmospheric 
distribution, related to a particular planet-star interac- 
tion scenario, may hide or compensate the loss signature 
during the observing time. Future HST (when repaired) 
FUV observation of the system during a full planetary 
orbit should help to disentangle the different processes 
in play. 
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